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ABSTRACT 


yA  class  of  ocean  acoustic  wave  propagation  problems  is  represented  by  a 
parabolic  equation  of  the  Schrodinger  type.  Using  conventional  explicit 
finite  difference  schemes,  e.g.,  the  Euler  scheme,  to  solve  the  parabolic  wave 
equation  is  unstable.  Thus,  important  advantages  of  explicit  schemes  are 
completely  missing.'  This  paper  presents  a  conditionally  stable  explicit 
scheme  by  introducing  an  extra  dissipative  term.  This  new  explicit  scheme  is 
then  applied  to  solve  the  ocean  acoustic  parabolic  wave  equation  fully 
utilizing  the  advantages  of  explicit  schemes.  The  theoretical  development, 
the  computational  aspects,  and  the  advantages  are  discussed.  Application  of 
the  scheme  to  a  realistic  ocean  acoustic  problem  is  included.  The  solution 
obtained  is  compared  with  the  unconditionally  stable  Crank-Nicolson  solution. 


INTRODUCTION 


A  parabolic  equation  of  tbe  Scbrodinger  type  arises  in  the  application  of 
ocean  acoustic  wave  propagation.  In  the  published  literatures,  this  ocean 
acoustic  parabolic  wave  equation  has  been  solved  by  three  different  methods: 
the  Split-step  Fourier  a1gorithm£^3,  the  numerical 
ordinary-differential-equation  raethodC^*^*^^,  and  an  implicit 
finite-difference  methodf5»6j.  The  implicit  finite-difference  method  is,  by 
far,  the  most  general  purpose,  stable  method  for  solving  the  parabolic  wave 
equation;  in  addition,  it  has  a  variety  of  useful  capabilities.  However,  when 
3-dimensional,  as  well  as  high  frequency  problems  arise,  one  needs  a  more 
effective  method  to  ease  the  requirement  of  memory  storage,  to  gain 
computational  speed,  and  to  be  implemented  easily  into  modern  pipe-line 
computers.  It  Is  known  that  explicit  schemes  have  these  advantages.  Way  back 
*  to  1978,  Lee  and  PapadakisL^j  analyzed  applicable  explicit  schemes  for  such 
application  and  found  that  the  explicit  scheme  such  as  the  Euler  scheme  was 
unstable'.  .'/It  is  the  main  result  of  this  paper  to  introduce  a  stable  explicit 
scheme,  newly  developed  by  Chan,  Lee,  and  Shenl^^3  for  solving  the  parabolic 
equation  of  the  Schrodinger  type.  This  new  explicit  scheme  Is  developed  to  be 
conditionally  stable  by  adding  a  dissipative  term,  and  possesses  most 
advantages  any  explicit  scheme  must  have.  Prior  to  the  discussion  of  the 
formulation  and  the  stability  of  the  new  explicit  scheme,  two  sections  are 
presented.  One  is  the  introduction  of  the  problem  background,  the  other  is 
the  discussion  of  the  solution  background.  An  application  to  a  realistic 
ocean  acoustic  problem  Is  presented  along  with  a  comparison  with  the  solution 
obtained  by  the  Crank-Nicolson  method  on  accuracy  and  speed. 
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PROBLEM  BACKGROUND 


In  order  to  simplify  compncated  ocean  acoustic  wave  propagation 
problems,  Tappmrt^^^  introduced  the  Parabolic  Equation  approximation  (PE)  to  v 

solve  the  elliptic  wave  equation  In  2  dimensions  by  a  parabolic  equation 
involving  the  depth  variable  z  and  the  range  variable  r.  This  PE 
approximation  produces  the  parabolic  wave  equation  in  the  form 

“r  ■  I  ^  “zz 

where  u(r,zj  is  the  pressure  fielo,  is  a  reference  wave  number,  and  the 
n(r,z)  is  the  inoex  of  refraction  which  is  equal  to  CQ/c(r,z)  (reference 
sound  speed/sound  speed).  Eq.  (1)  is  •  .yaroed  as  the  standaro  parabolic  wave 
equation  first  introouced  by  TappertC^].  , 

The  parabolic  wave  equation  (1)  can  be  obtained  in  a  couple  of  different 
ways.  We  give  an  outline  of  our  development  below. 

The  elliptic  wave  equation  in  cylindrical  coordinates  takes  the  form 


ar  az^^ 


where  b(r,z)  is  the  acoustic  pressure  field. 


(2) 


Expressing 


where  v(r)  has  strong  dependence  on  the  range  variable  r  while  u(r,z)  is 
weakly  dependent  on  r.  We  shall  derive  that  the  u(r,2)  here  satisfied  Eq.  (1). 

Substituting  expression  (3)  into  Eq.  (2),  we  obtain 


[Vr  *  7  *r]  “  *  [“rr  *(7*7  *r)  “r  ’  “j2  ’  "0  ^  •  » 


(4) 
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Setting  the  terms  in  the  first  [  ]  of  Eq.  (4)  equal  to  *-k  v*  and  setting  the 
terms  In  the  second  [  ]  of  Eq.  (4)  to  "k^u,"  we  find  ® 


1  ? 

V  ♦  i  V  +  k  V  ■  0  .  IS) 

rr  r  r  0  “  *  '  ^ 

and 

“rr  *  (7  *  7  “r)  “r  *  “zz  *  ‘5("^('-.z)-l)u  -  0  (6) 

Consioering  only  the  outgoing  wave  In  the  range  direction,  it  is  easily 
seen  that  the  solution  of  Eq.  (5)  Is  the  zeroth  order  Hankel  function  of  the 
first  kino.  If  we  apply  the  far-fielo  approximation,  k^r  »  1,  to  the 
argument  of  the  Hankel  function,  we  obtain 


v(r)  - 


i(V-t) 


(7) 


Expression  (7)  Is  used  to  simplify  the  coefficient  (p 7  v^)  of  Eq.  (6) 
which  becomes 

“rr  *  21Vr  “zz  u  .  0  .  (8) 
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Making  use  of  the  property  of  n(r,z}  being  slowly  varying  In  r,  and 
neglecting  the  scattering  In  both  directions.  Eg.  (8)  can  be  rewritten  In  an 
operator  form 


Since  a  one-way  outgoing  wave  Is  considered,  we  can  deal  with  the 
solution  of  the  out-going  wave  equation 


Applying  a  low  order  approximation  to  the  square-root  operator  In  Eq.  (10),  we 
obtain 


2 

wil  ♦  ■^(n^-1)  *  ^  •  (11) 

kjj  3Z 

Substituting  (11)  Into  Eq.  (10) and  simplifying,  the  following  equation  Is 
obtained 


A  ♦  (n2-l) 

'  kf  az2 


■  7  2^  “zz 


which  is  exactly  the  same  as  Eq.  (1). 


It  Is  understood  that  If  a  solution  is  sought  for  the  elliptic  wave 
equation,  Eq.  (2),  one  must  solve  a  purely  boundary  value  problem  in  a  region 
as  described  below  in  Figure  1. 

where  is  the  initial  boundary  condition, 

#('‘(2^)  IS  the  surface  boundary  condition, 

0('‘f2g)  Is  the  bottom  bounoary  condition,  ano 
is  the  wall  boundary  condition. 

On  the  other  hand,  if  the  elliptic  problem  was  solved  by  the  parabolic 
equation,  Eq.  (1),  we  deal  with  a  solution  of  initial  boundary  value  problem  - 
as  shown  by  Figure  2  where  we  need  only  to  know  ulr^.z)  the  initial 
conoition,  u(r,2^)  the  surface  boundary  condition,  and  uCr.zg)  the  bottom 
boundary  condition  while  the  wall  boundary  can  be  ignored  completely. 


Figure  1.  Region  of  Elliptic  Figure  2.  Region  of  Parabolic 

Problems  Problems 

Immediately,  we  see  the  advantage  that  the  parabolic  problem  requires  one  less 
^  boundary  condition  on  the  wall  which  Is  usually  difficult  to  specify.  In 
addition,  the  parabolic  problem  can  be  solved  by  a  marching  process  stepping 
forward  in  range.  It  would  be  more  advantageous  if  the  parabolic  problem  can 
be  solved  very  efficiently.  This  is  the  major  result  we  discuss  in  this  paper. 


SOLUTION  BACKGROUND 


Three  different  methods  exist  to  solve  the  parabolic  wave  equation,  Eq. 
(1).  These  methods  are:  Split-step  algorithmt^^,  numerical 
ordinary-olfferentlal-equatlon  inethodsL^ »3 .4] ^  ^  implicit  Finite 

Difference  (IFO)  methodt®*®^.  Ths  split-step  Fourier  algorithm  Is  effective 
for  deep  water  problems  where  the  pressure  field  vanishes  at  both  surface  and 
bottom  boundaries.  The  numerical  orolnary-dlfferential-equatlon  solution 
applies  the  Generalized  Adams  methoqt^j,  recently  much  Improveo  by  Lee, 
Jackson,  ano  PreiserC4].  The  method  is  effective  and  general  purpose  and  It 
can  handle  arbitrary  boundary  conditions.  Most  frequently  used  In  the 
application  of  ocean  acoustic  propagation  for  medium  to  low  frequency  problems 
Is  the  implicit  finite  difference  (IFD)^®]  solution.  This  implicit  finite 
difference  (IFO)  solution  Is  not  only  general  purpose,  unconditionally  stable, 
but  also  has  wider  angle  capabilities  than  other  existing  methods.  Trying  to 
take  advantage  of  requiring  less  storage,  ano  the  ease  of  implementation  into 
pipe-line  computers,  Lee  and  Papaoaklsi-^^  maoe  an  attempt  to  apply  explicit 
schemes  such  as  the  Euler  scheme  to  solve  the  parabolic  wave  equation.  Their 
analyses  showeo  that  using  the  Euler  explicit  scheme  to  solve  Eq.  (1)*  the 
scheme  is  unstable.  Then,  Lee,  Botseas,  and  Papadakis  applied  the  IFD,  which 
uses  the  Crank-Nicolson  scheme  to  solve  Eq.  (I)  because  of  Its  favorable 
unconditional  stability.  Even  though  the  stability  is  maintained  by  the  IFO 
scheme,  the  desirable  advantage  of  explicit  schemes  are  completely  missing. 
Recently,  Chan,  Lee,  and  ShenL^i  developed  a  group  of  explicit  schemes  for 
solving  the  equations  of  the  Schrodinger  type.  A  new  explicit  scheme  to  be 
introouceo  in  the  next  section  is  a  member  of  that  group  which  is  formulated 
to  be  conditionally  stable  while  other  important  advantages  of  the  explicit 
scheme  are  all  retained. 
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A  CONOITIONALLY  STABLE  EXPLICIT  SCHEME 


Compared  to  implicit  scDemes,  explicit  schemes  are  generally  easier  to 
implement  and  demand  less  storage.  These  advantages  are  especially  pronounced 
for  multi-dimensional  problems.  Moreover,  another  important  advantage  is  that 
the  explicit  scheme  is  often  easily  vectorized  on  the  many  pipeline-oriented 
computers  availaole  today  such  as  Crays,  Cyber  205  and  FPS164.  Talcing  these 
advantages  our  method  of  attacic  is  to  introduce  an  appropriate  dissipative 
term  to  derive  stable  explicit  schemes. 

We  begin  by  considering  the  parabolic  equation  of  Schrodinger  type  in  the 
general  form  below 

Up  •  ia(r,2)  U22  +  b(r,z)  u^  ♦  c(r,2)  u  ♦  f(r,2)  (12) 

Here  the  function  a(r,z}  is  a  real-value  function,  however,  functions  b(r,z), 
c(r,z),  and  f(r,z)  may  be  complex-valued.  The  effects  of  adding  different 
dissipative  terms  have  been  discussed  in  Ref.  7.  In  this  paper  we  discuss  one 
particular  dissipative  term  which  leads  to  the  least  restrictive  stability 
condition.  We  use  k  ano  h  to  denote  the  range  and  depth  Increments 
respectively,  u^  means  u(r",  z^). 

The  simplest  explicit  scheme  for  a  simple  Schrodinger  equation, 

“r  -  1u„  .  (13  ) 


is  the  Euler  scheme. 
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(14) 


with  an  initla)  truncation  error  0(k»h2}.  Scheme  (14)  Is  known  to  be 
unstable  We  Introduce  a  dissipative  term  R  to  Eq.  (13)  to  obtain 


V  -  luzz  *  ^  <15) 

where  the  R  In  this  paper  Is 

R  -  (a-^lB)  h2  U^2ZZ  (16) 

where  a  and  a  are  real  scalars.  With  the  dissipative  term  R,  the 
corresponding  explicit  scneme  for  Eq.  (13)  is: 


.  *  "j-i^  . 


-  4u . 


♦  6u' 


■  “i-l- 


(17) 


ano  the  least  restrictive  stability  constraint  Is 


^  .  1 


which  Is  obtained  when  a 


1 

-  V 


a  •  X 


(18) 


To  show  that  scheme  (17)  Is  stable  when  It  Is  used  to  solve  Eq.  (13),  we  quote 
a  theorem  below  which  has  been  proved  In  Ref.  7. 
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(19) 


THEOREM:  The  scheme  (17)  is  stable  if  and  only  if  o  <  0  and 

min /-2a,  - y — - 

h^  \  16a^  +  (4e-i)V 

Tne  definition  of  stability  used  here  is  the  notion  of  practical  stability  as 
discussed  by  RichtmeyeWlorton  Lflj  and  Chan  [9]  which  requires  the  discrete 
solution  to  have  a  nongrowing  norm. 

Now  we  extend  scheme  (17)  to  solve  the  more  general  equation,  Eq.  (12). 
It  is  easily  seen  that  the  extenoeo  stable  explicit  scheme  takes  the 
expression 


Under  a  slightly  weaker  stability  definition,  the  stability  of  this  extended 
scheme  is  'Iso  given  by  (18). 

He  have,  thus,  introouceo  a  conditionally  stable  explicit  scheme  with 
least  restrictive  stability  condition.  We,  now,  proceed  to  show  its 
advantages.  We  begin  by  applying  this  scheme  to  a  real  problem  ano  comparing 
its  results  with  the  results  produced  by  the  Crank-Nicolson  scheme  in  the 
section  to  follow. 


AN  APPLICATION 


We  apply  the  explicit  scneme  (20)  to  solve  a  wave  propagation  proolem  in 
the  region  of  the  Mediterranean  Sea.  The  representative  wave  equation  Is  Eg. 
(1).  This  is  a  propagation  problem  under  a  shallow  water  environment  where 
the  water  depth  Is  100  m.  Under  such  environment*  an  isovelocity  sound  speed 
Is  considered.  The  bottom  is  characterized  by  a  slightly  different  sound 
speed*  and  a  different  water  density  Is  specified;  In  addition*  bottom 
attenuation  Is  required  due  to  the  bottom  energy  loss. 

Both  the  source  and  the  receiver  are  placed  at  the  same  depth  In  the 
ffliodle  of  the  region.  The  source  propagates  at  a  low  frequency.  We  predict 
propagation  loss  up  to  the  range  of  25  kilometers.  This  problem  has  been 
solved  by  three  different  methods  -  a  normal  mode  method  (SNAP  [11 j).  a 
split-step  Fourier  algorithm  (PAREQ  [11]),  and  the  implicit  finite  difference 
methoo  (IFD  L6j).  The  solutions  produced  by  the  above  three  methoos  agree 
exactlyi-^J.  The  SNAP  normal  mode  solution  was  used  as  a  benchmark  reference 
solution  for  comparison.  The  PAREQ  and  IFD  solutions  are  all  performed  in  a 
marching  process*  In  oroer  to  satisfy  the  Fourier  requirements*  the  PAREQ 
used  512  points  as  a  depth  partition,  then  marcheo  with  a  range  step  of  1/2 
meter.  The  IFD  useo  exactly  the  same  step  sizes  in  oroer  to  make  a 
point-to-point  comparison.  Results  are  all  satisfactorily  accurate,  however* 
the  IFD  Is  two-thirds  faster  than  the  PAREQ.  For  the  discussion  of  the 
application  of  this  new  stable  explicit  scheme*  we  choose  to  compare  with  the 
f Inite-olfference  (IFD)  solution  In  terms  of  accuracy*  speed*  and  Its 
implementation  effort. 


In  solving  this  problem,  the  following  Input  parameters  were  used: 


source  depth 

-  50  m 

source  frequency 

■  25  Hz 

bottom  depth 

-  100  m 

receiver  depth 

-  50  m 

souno  speed  (water) 

>  1500  m/s 

sound  speed  (bottom) 

.  1550  m/s 

density  (water  column) 

«  1  g/cm^ 

density  (at  bottom) 

■  1.2  g/cm3 

attenuation 

-  1  dB/wavelength 

maximum  range 

■  40  km 

Initial  range 

-  zero  m 

(  The  result  obtained  Is  Propagation  Loss  (PL)  ano  Is  measured  In 

"decibels,"  abbreviated  "dB."  Propagation  Loss  In  underwater  acoustics 

« 

terminology  quantitatively  describes  the  weakening  of  sound  as  It  travels 
through  the  sea.  The  conventional  formula  for  the  propagation  loss  Is  PL  ■ 
-20  1og2g(|u|).  A  graph  of  PL  vs  Range  Is  plotted,  as  shown  In  Figure  3 
where  both  IFO  and  explicit  numerical  results  are  Identical  to  one  another. 
The  explicit  scheme  produced  satisfactorily  results,  however,  the  computation 
speed  of  the  explicit  scheme  is  a  little  faster  than  the  IFD  method  (114  min 
IFO  vs  90  min  explicit  scheme).  The  largest  allowable  range  step  size 
required  by  the  explicit  scheme  for  stability  was  ~  0.5  m,  which  agrees  with 
the  theoretical  estimate  (18). 
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RANOSillW) 


Figure  3.  Solutions  Comparison 


CONCLUSIONS 

A  conoitlonall^  stable  explicit  finite  difference  scheme  which  can  be 
applied  to  shallow  water,  low  frequency,  long  range  wave  propagation  problems 
In  the  ocean  has  been  Introouceo.  Comparison  showeo  the  advantage  of 
computational  speed  over  the  Crank-NIcolson  solution  because  the  explicit 
scheme  required  no  solution  of  a  system  of  equations,  but  scalar  operations. 
This  advantage  automatically  Indicates  the  relaxation  of  memory  storage.  It 
also  Indicates  an  easy  vectonzatlon  of  the  explicit  scheme  on 
pipeline-oriented  computers. 

The  scheme  reported  In  this  article  Is  a  member  of  a  group  of  stable 
explicit  schemes.  These  schemes  are  primarily  developed  to  solve  the  equation 
of  the  Schrooinger  type  in  a  general  nature.  It  is  expected  that  this  group 
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of  explicit  schemes  will  have  equally  efficient  applications  of  problems  in 
other  fields  such  as  plasma  physics,  quantum  mechanics,  and  seismology  due  to 
their  basic  advantages  -  fast,  requiring  less  memory  storage,  easy  to 
implement,  and  easy  to  vectorize  on  pipeline-orienteo  computers. 
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